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Quantum transport through disordered structures is inhibited by (Anderson) localization effects. The disorder 
can be either topological as in random networks or energetical as in the original Anderson model. In both cases 
the eigenstates of the Hamiltonian associated with the network become localized. We show how to overcome 
localization by network multiplexing. Here, multiple layers of random networks with the same number of nodes 
are stacked in such a way that in the perpendicular directions regular one-dimensional networks are formed. 
Depending on the ratio of the coupling within the layer and perpendicular to it, transport gets either enhanced 
or diminished. In particular, if the couplings are of the same order, transport gets enhanced and localization 
effects can be overcome. We exemplify our results by two examples: multiplexes of random networks and of 
one-dimensional Anderson models. 

PACS numbers: 05.60.Gg, 05.60.Cd, 71.35.-y 


Introduction - On a coarse grained scale, many physical, 
chemical, biological, or sociological systems can be modelled 
by networks of connected nodes, where every node represents 
one of the partners (atoms, molecules, people, etc.) interact¬ 
ing with each other via the bonds. Recently, there has been 
active research on networks of networks, i.e., on possibly dif¬ 
ferent networks which are interdependent EKD- A particular 
example are so-called multiplex networks, which are layers 
of interdependent networks BHD, see also Fig. |T] One of 
the interesting aspects is the (incoherent/diffusive) dynamics 
of excitations within these networks Q, especially when the 
networks within a layer are non-regular, e.g., random or of 
small-world type (8). 

When one is concerned with the coherent quantum dynam¬ 
ics of an excitation on a potential landscape given by the net¬ 
work, the dynamical properties can be vastly different El, 
while the static properties of networks can be considered to 
be the same also for quantum mechanical problems. There 
are several ways to quantify the (global/averaged) efficiency 
of the coherent and of the incoherent dynamics, one of which 
being the so-called first passage time, i.e., the time it takes for 
an excitation to first reach a given node of the network. If this 
node is equipped with a decay channel such that the excitation 
can leave the network, the time integral of the survival prob¬ 
ability of the excitation defines the mean first passage time 
(MFPT) liTOl . This definition can be applied for both, inco¬ 
herent and coherent processes, where the latter needs some 
adjustment, see below. For diffusive processes such as ran¬ 
dom walks, MFPTs have been calculated, e.g., for polymer 
reaction kinetics ED or confined systems lfl2l . see Ifl3l for 
an overview of different applications. 

For the quantum case, non-regular networks can lead to lo¬ 
calization effects, e.g., if the network (in each layer) can be 
described by a random matrix ED or by a one-dimensional 
Anderson model E3- These localizations prevent coher¬ 
ent transport, consequently leading to very large first passage 
times. As we are going to show, this effect can be overcome 
by considering multiplex networks of non-regular intra-layer 
networks which are regularly interconnected, see Fig. [T] 


Concept - We take the idea of multiplex networks to the 
quantum domain. To each node k of the network one asso¬ 
ciates a basis vector/state | k), all of which form an orthonor¬ 
mal basis of the available state space, i.e., the Hilbert space 
for quantum mechanical problems El- The connectivity of 
the network is captured by the adjacency matrix, which has 
only non-zero elements if two nodes are connected by a sin¬ 
gle bond, implying that the diagonal elements are zero. The 
non-zero elements can be all equal, as, e.g., in random net¬ 
works for which the distribution of the non-zero elements is 
random lf8l H6l . hence, becoming a random matrix fl4l . It is 
also common to have non-zero elements not being equal, but 
rather randomly distributed in value with a regular distribu¬ 
tion in the adjacency matrix, as, e.g., in the Anderson model 
El- In any case, we will identify the adjacency matrix with 
the (discrete) Hamiltonian of the system, see below. 


decay 



FIG. 1. (Color online) Sketch of a multiplex network with (random) 
networks of N nodes in each of the M layers. 

We are further interested in the excitation transport through 
such networks. By allowing for localized decay of the prob¬ 
ability at specified nodes, we are able to monitor the global 
decay of the survival probability of the excitation which is a 
measure for the transport efficiency El M El: If the decay 
is fast, transport from an initial node to the node from which 
decay is possible is efficient and vice versa. 

For each layer to (to = 1,..., M) of the multiplex net- 
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work, we label the states corresponding to the nodes by 
|l m ),..., | N m ). When having (localized) decay, the prob¬ 
lem can be treated by employing an effective non-Hermitian 
intra-layer Hamiltonian m which reads 

H ™= E v LkJjm)(k m \~ir m \N m )(N m \ (1) 

Umk m ) 

where vj k is the coupling strength between pairs (j m k m ) 
(jm ) k m = 1,..., N m ) of nodes which are coupled by a sin¬ 
gle bond and where r m is the decay rate from node \N m ) (in 
the following we will assume T m = F for all in). There¬ 
fore, the intra-layer dynamics follows for each layer m from 
the Liouville-von Neumann equation for the density operator 
Pm (f) for a single layer m, Pm (f) — Pm(()] 11811191 . 

Note, that the right-hand side of this equation has a “non co¬ 
herent” part arising from the decay terms in Hm. 

Now, every node j rn in layer m is coupled to a single corre¬ 
sponding node j m +i in the adjacent layer m+ 1 with coupling 
strength V^ m+1 . This creates chains of nodes in the direction 
perpendicular to the layers. Then, the inter-layer Hamiltonian 
is given by 


where A/o is the number of vanishing values of the 7 „’s. 
Therefore, and in analogy to the classical case flOl . we now 
define the quantum analog of the global mean first passage 
time (gMFPT) as 


rs /dt[n(t)-nj = T £ 2-. (5) 

0 {«l7n^0} 

Clearly, Eq. ([5]) depends on Wo as well as on N and on M. In 
the following we will keep N constant while varying M. 

Since the intra-layer networks we consider are random net¬ 
works we will further analyze the ensemble averages over 
R realizations of the multiplex network, namely (11 ^« = 
Ef=i nS and ( t) r = where the superscript 

[r] denotes the rth realization. In order to obtain analyti¬ 
cal approximations/bounds, we calculate the quantum analog 
of the gMFPT from the ensemble averaged 7 n ’s, (7 h )r = 
Eti ln\ i.e. we define 


T 


\ 
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E 


1 

2(7 u)r 


( 6 ) 


M—1 

ffj = E! En,m+lbm)Om*l|- (2) 

m= 1 

The full dynamics for the density operator p(t ) of the whole 
multiplex network is then described by p(t) = —i\H. p(t)\ 

with H = J2m =1 + SjLi H ) In total there are M 

nodes from which an excitation can decay, with corresponding 
states |W m ) with m = 1,..., M. 

In order to study the transport through these networks, we 
assume initial conditions which are localized at any node of 
the whole multiplex network such that the transition proba¬ 
bility from node k to node j reads nj k (t) = (j\p(t)\j) with 
p(0) = |fc)(fc|. The problem can be properly solved by diago- 
nalization of H. In most cases, the eigenvalues E n = e rl +T/ rl 
(n = 1 ,.. •, NM) are complex with e n £ K. and 7 „ £ Rq\ 
By then averaging over all possible initial nodes k and all pos¬ 
sible final nodes j we obtain a global picture of the decay pro¬ 
cess, i.e., of the transport efficiency. Thus, the mean survival 
probability is given by naira: 

1 NM NM 

n W ~nm^ njk{t) = NM £ ex P(-27nf). (3) 

k,j =1 n—1 

In the limit when F < 1, the imaginary parts are related 
to the eigenstates of the network without decay, namely 
7 n = r E m |(-^m|$r l 0) )| 2 , where |$ 1 0) ) are the eigenstate of 
if(°) = Hm 0 ' 1 + Hj~, where Hm °' ) is the hermitian 

part of Hm 0 . 

General Results - Now, since not necessarily all values of 
the 7 n ’s need to be larger than zero, the mean survival proba¬ 
bility can decay to a finite value given by 

.^ n <‘> sn ~ = M? E ' 4 ' (4) 

{n| 7 „= 0 } 


Jensen’s inequality guarantees that r < (t)r, see paragraph 
12.41 of Il20l. As will be shown in what follows, the properties 
of T and (r)n, as well as of (No)r crucially depend on the 
number of layers and on the ratio V" /V L . For the two limits 
V 7> V I and V /X <C perturbation theory yields that T 
and also (■ t)r decrease with increasing M: 

(i) For V 1 - F b the total networks consists of N chain¬ 

like networks, each of length M, which are only weakly cou¬ 
pled within each layer. This implies that approximately only 
those states |$(°)) belonging to the chain which contains the 
nodes with decay have non vanishing overlap with the states 
| N m ). Therefore, for /V 1 - 0 there will be only M 

values 7 „ ^ 0 entering r and consequently also only M 
averaged values ( 7 h )r ^ 0 entering t. Further, one finds 
Ao ~ (N—1)M. The eigenstates of each chain are Bloch-like 
states, i.e., they are states of the form (1/y/M) E m c m\N rn ), 
where c m is a complex number lying on the unit circle ED. 
This then leads to r ~ (1/2JVT) E m l/|cm| 2 , with c m = 1 
for Bloch-like states, such that r ~ 1 / (2FAT) becomes inde¬ 
pendent of M. This result also holds in the ensemble averages, 
both, for r and for (t)r. 

(ii) In the other limit where V - 1 <C V^, the total multi¬ 
plex network consists of very weakly coupled layers. Thus, 
the dynamics is mainly dominated by the dynamics within a 
single layer and in addition by weak tranport perpendicular to 
this layer. In this case, however, we have to consider the en¬ 
semble averages since the details of each realization of a ran¬ 
dom network can vary (strongly). The (average) eigenstates 
of random networks are exponentially localized (H, such that 
a significant overlap of these states with | N m ), occurs only 
for a number of order 0(1) of these states. When assuming 
for simplicity that all the other states have zero overlap, the 
sum in r will only run over those number of states with over¬ 
lap which for the whole multiplex network is of order 0{M), 
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again implying that TVq ~ (N — 1 )M. Consequently, also 
here r is independent of M. As it turns out, the same holds 
for the ensemble average (' t)r■ 

Having, in the two limits (i) and (ii), for an increasing num¬ 
ber of layers, M, a linear increase in the number of vanishing 
values of the 7 n ’s and constant values for the quantum analogs 
of the gMFPTs does not imply that this is true for all values 
of V''/V 1 . As our calculations show, one finds a decrease of 
all considered quantities with M when V~ l « V", i.e., when 
no perturbation limit can be applied. 

The two limits (i) and (ii) show that the minimal num¬ 
ber of values 7 „ = 0 is given by J\f Q ~ (N — 1 )M. For 
fixed M, any change in the couplings leading away from 
these limits will lead to an increase of this number. This is 
best seen on the basis of the eigenstates |4>^): A change 
in the ratio V''/V 1 - away from the limits (i) and (ii) re¬ 
sults in eigenstates which are distributed over larger areas of 
the multiplex network. Consequently, the number of values 
7 „ = T |(iV m |$i°^)| 2 = 0 will decrease and more prob¬ 
ability can leak out of the system, leading to a lowering of 
11^. Now, increasing M will result in a further decrease of 
Afo and of n^. 

Similar arguments hold for r: Since the eigenstates be¬ 
come “broader” with in-/decreasing the ratio V'' /V L from 
the limit (i)/(ii), the overlap of most eigenstates with \N m ) 
becomes larger, leading to a larger value of the corresponding 
7 This happens on the cost of the overlap of other eigen¬ 
states with \N m ), since Y2 n l(-^ r m| < £’ri 0 ^)| 2 = 1 holds. Sup¬ 
pose that there are only two states with overlap such that, 
say, dj = (Ay„|<Iv 0 ' l )| 2 with j = 1,2, for which ai > 02 
and 01+02 = 1. After a change in we assume 

to have new states with di = ai — e and a 2 = 02 — e, 
where 1 e > 0, such that also 01 + 02 = 1. For r 
we need 1/di + 1 / 02 , which after expanding in e becomes 
l/oi + 1/02 - e(l - 2o 2 )/[a|(l - a 2 )j < 1/ai + l/a 2 . As a 
result, these change in the values of the 7 „’s leads to a smaller 
value of r. Also here, increasing M leads to a further decrease 
of r. 

Examples - We will exemplify the above results by two 
variants of multiplex networks. In order to minimize the nu¬ 
merical effort but still obtain reasonable results, we take for 
both examples N = 10, a trapping rate of T = 0.001 • V", 
and the number of realizations R = 200. 

First, we consider multiplex networks where each layer 
consists of a random network of given size N, as sketched 
in Fig. [] The maximal number of bonds in each network 
is N(N — l)/2, out of which we take a fraction of 1/2 to 
be placed at random between the bonds, such that maximally 
one bond connects two nodes. We assume the same coupling 
strength V" = 1 for any pair of nodes which is connected by 
a single bond. In the second example, we take the network in 
each layer to be a one-dimensional Anderson model with off- 
diagonal disorder, again with N = 10. Here the disordered 
couplings vj k = + Aj mt k rn are drawn from a normal 

distribution centered around the value V" =1. In both cases 


V ± will be given in units of V". 

Since all quantities we consider depend on the values of 
the 7 „’s, in particular on those values which are either exactly 
zero or strictly larger than zero, it makes it numerically dif¬ 
ficult to determine whether a value is truly zero or not. This 
can cause large deviations in, say, the values of r, see Eq. 0- 
Therefore, we use a cut-off below which we assume the 7 n ’s 
to be zero, i.e., we neglect long-time tails of n(f) and also 
slightly raise the value of 11^. We have checked different 
values of the cut-off and found a value of 10 h a reasonable 
choice, because the fluctuations in the considered quantities 
were rather small, see the numerical results below. We note 
that choosing too large a cut-off could deminish the differ¬ 
ences discussed above. 


(a) (n oo ) H /max[(n oo ) iJ ] (b)(r) B / max[(r)i?] (c) r/max[r] 
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FIG. 2. (Color online) Multiplex networks with each layer being 
a random network of size N = 10 with a fraction of 1/2 of all 
possible bonds present, (a)-(c) show contour plots of (noo).R, {t)r, 
and t, respectively, for different ratios of the couplings H^/V r± and 
different numbers of layers, M. Note that all plots are scaled by the 
respective maximal values. There is a clear minimum for all three 
quantities at value H^/H x fa 1. The lower panels illusbate the in¬ 
dependence (d) of (noo)_R for different values of H^/H x and (e) of 
{t)r and t for H^/H x = 1. Especially for (t)r there is a drop by 
about one order of magnitude from M = ltoM = 10. 


For multiplexes of random networks, Figs. |2|a)-(b) show 
(noo)ii. (t)r, and r, respectively, for different ratios [V L 
and different M and rescaled to their respective maximal val¬ 
ues. The first thing to notice is that in all three cases there 
is a pronounced minimum for ratios V'H / V 1 fa 1, which 
decreases with increasing M, thus corroborating the general 
statements given above. This is further illustrated in Fig.[2jd) 
where the M-dependence of (11^) for different values of 
V' 11 / V 1 is shown and in Fig.^e) where the M-dependence of 
(r) r and r for / V ± = 1 is displayed. In particular, (t)r 
shows a drop by about one order of magnitude from M = 1 
to M = 10. This indicates that one can overcome localization 
effects by a suitable choice of multiplexing several (random) 
networks. Moreover, the results also confirm that in the two 
limits (i) and (ii) all three quantitites do not change with M 
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for M > 2. For all ratios V''/V 1 - , there is a change in all 
quantitites when going from M = 1 to M = 2, see Fig[2fd) 
and (e). In addition, we find the trivial result that for M = 1 
(only one layer, i.e., no multiplex) the three quantities do not 
change with the ratio IZlI /V^, see Fig§d). 


(a) (n oo ) B /max[(n oo ) fl ] (b)(r) B / max[(r)i;] (c) r/max[r] 
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FIG. 3. (Color online) Multiplex networks with each layer being an 
Anderson model of N = 10 nodes with off-diagonal disorder drawn 
from a normal distribution with width a = 0.5V^. See Fig. [ 2 ] for 
description of the different panels. 


For multiplexes of Anderson models, Fig.[3]shows the same 
quantities as Fig. [2] for slightly different parameter ranges of 
V''/V 1 - and M. Given similar features of the eigenstates 
|$(°)) of random networks and of the Anderson model, we 
find similar results as for the multiplexes of random networks, 
see Figs. [3ja)-(c). Most importantly, also here there is a pro¬ 
nounced minimum at values / V ± ~ 1, where (IIoo^ de¬ 
creases by about two orders of magnitude from M = 1 to 
M = 6, see Fig. [3jd). The quantum analog of the gMFPT, 
( t)r , decreases in the same interval of M by about one order 
of magnitude. Again, suitable multiplexing of several Ander¬ 
son model overcomes localization and the overall transport 
through the whole network becomes (much) more efficient. 

Conclusions - We have shown that the coherent quantum 
transport on non-regular networks which show localization 
can be enhanced arranging these networks in layers which 
are regularly interconnected. The quantum analog of the 
global mean first passage time allows to define a global mea¬ 
sure for the transport efficiency, which depends on the ratio 
between the intra- and the inter-network coupling strengths, 
/V 1 - and the number of layers, M. The ratio / V ± de¬ 

termines how large the enhancement will be: While in the 
limits /V 1 - < 1 and IZlI /V 1 - 1 there is only very little 

enhancement, one can achieve a significant enhancemnt for 
V''/V 1 - ss 1. Our general analytical results are corroborated 
by numerical calculations for two examples, multiplex net¬ 
works of intra-layer random networks and of intra-layer one¬ 


dimensional Anderson models. We believe that these model 
systems can also be realized experimentally, say, using wave¬ 
guide arrays for which, e.g., Anderson localization has been 
realized I22ti24l . One can also envision that multiplex net¬ 
works could be used to improve signal transfer. 
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